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ABSTRACT 

Fluorescent nucleic acid base analogues are im- 
portant spectroscopic tools for understanding 
local structure and dynamics of DNA and RNA. We 
studied the orientations and magnitudes of the 
electric dipole transition moments (EDTMs) of 
6-methyl isoxanthopterin (6-MI), a fluorescent 
analogue of guanine that has been particularly 
useful in biological studies. Using a combination of 
absorption spectroscopy, linear dichroism (LD) and 
quantum chemical calculations, we identified six 
electronic transitions that occur within the 25 000- 
50000 cm -1 spectral range. Our results indicate that 
the two experimentally observed lowest-energy 
transitions, which occur at 29687 cm -1 (337 nm) 
and 34596 cm -1 (289 nm), are each polarized within 
the plane of the 6-MI base. A third in-plane 
polarized transition is experimentally observed at 
47547 cm -1 (210 nm). The theoretically predicted 
orientation of the lowest-energy transition moment 
agrees well with experiment. Based on these 
results, we constructed an exciton model to 
describe the absorption spectra of a 6-MI dinucleo- 
tide-substituted double-stranded DNA construct. 
This model is in good agreement with the experi- 
mental data. The orientations and intensities of the 
low-energy electronic transitions of 6-MI reported 
here should be useful for studying local conform- 
ations of DNA and RNA in biologically important 
complexes. 



INTRODUCTION 

Experiments that reveal the local conformations of nucleic 
acids and conformational transitions along a reaction 
pathway are needed to better understand the mechanisms 
of the various proteins and protein complexes that 
comprise the 'macromolecular machines' that manipulate 
genes and other elements of the DNA genome. The use of 
fluorescent nucleic acid analogue probes, in conjunction 
with well-chosen spectroscopic measurements, can help to 
reveal important information about local conformations 
at chosen sites within nucleic acid constructs, and thus the 
roles played by these specific structures in key steps of 
biologically important protein-nucleic acid reactions. 

6-methyl isoxanthopterin (6-MI) is a fluorescent 
analogue of the natural nucleic acid base guanine. 6-MI 
absorbs light in the long wave ultraviolet (UV) regime at 
lower energies (X ex ,max ^ 350 nm) than do the canonical 
nucleic acid bases, and (unlike the natural bases) 
fluoresces with relatively high quantum yield 
{Km,ma\ ^ 430 nm, Q » 0.7 for the nucleoside in aqueous 
buffer) (1-9). The ability of 6-MI to form Watson-Crick 
base pairs with cytosine when it replaces guanine within a 
DNA or RNA strand, in addition to its low-energy 
spectral properties, makes it a useful spectroscopic 
site-selective probe of local base conformation and 
dynamics (3,10,11). As is seen for other fluorescent 
nucleic acid base analogues such as 2-aminopurine, 
pyrrolo-cytosine and tricyclic cytosine, the spectral 
properties of 6-MI can provide information about the 
local conformation and dynamics of nucleic acids 
(3,9,12-15). The extent of fluorescence quenching in 
6-MI appears to be a measure of the degree of stacking 
between the site-specifically inserted base analogue probe 
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and its flanking bases on the same oligonucleotide strand. 
It has recently also been demonstrated that the substitu- 
tion of two sequential 6-MIs for adjacent guanines within 
DNA constructs results in exciton coupling between these 
'dimer' probes characterized by distinguishing features in 
the near-UV absorption and circular dichroism (CD) 
spectra (3,11). 

Exciton coupling is a spectroscopic phenomenon in 
which electronically interacting chromophores are 
excited collectively, with the electronic transition 
energies of the interacting components being split 
relative to those of the isolated molecules (16). The tran- 
sition energies and intensities of the electronically coupled 
dimer are sensitive to the relative orientation and separ- 
ation of the monomer subunits. Spectroscopic methods 
that are sensitive to exciton coupling, such as CD and 
2D fluorescence spectroscopy (2D FS) (17-19), can 
detect changes in the conformation of an exciton-coupled 
dimer. Thus by exploiting the electronic properties of 
6-MI dimer-substituted nucleic acids, local base-stacking 
conformations in various DNA (and RNA) constructs can 
be studied under physiological conditions (11). 

To better understand the spectroscopic properties of 
6-MI-substituted nucleic acids, it is important to charac- 
terize the electric dipole transition moments (EDTMs) of 
the isolated probe molecule. Interactions between the 
6-MI chromophore and its local environment perturb 
the electronic excitation and relaxation pathways, which 
is the origin of the probe's sensitivity to local macromol- 
ecular conformation. Knowledge of the directions and 
magnitudes of the EDTMs can be used to model the 
effects of nucleic acid structure on spectroscopic observ- 
ables such as transition strengths, energies and optical 
activity. Of particular importance for biochemical 
studies is the ability to model the electronically coupled 
states of two adjacent 6-MI chromophores substituted at 
sites within an oligonucleotide construct that are function- 
ally important in protein-nucleic acid complexes. Exciton 
coupling models have been used to elucidate local probe 
conformation in such constructs using CD spectroscopy 
(11,17). A combined analysis of linear spectroscopic meas- 
urements and 2D FS can provide more detailed informa- 
tion about the conformations of exciton-coupled 
chromophores than CD. For example, the conformation 
of a dimer of Mg-tetraphenylporphyrin suspended in a 
phospholipid bilayer membrane was recently determined 
using this approach (17,18). 

In the present work, we describe experimental and the- 
oretical investigations of the EDTMs of 6-MI. Although 
the photo-physical properties of 6-MI such as fluorescence 
parameters (1,7-9,20,21), absorption spectra (7,8,20,21), 
solvatochromism (8,20) and Stark spectra (8) have been 
the focus of previous experimental and theoretical studies, 
the directions and magnitudes of the EDTMs within the 
molecular frame have not been investigated experimen- 
tally. To determine transition polarizations and oscillator 
strengths, we performed UV absorbance and linear 
dichroism (LD) measurements on 6-MI nucleosides 
embedded in stretched poly(vinylalcohol) (PVA) films 
using an established approach developed by others 
(19,22-24). Comparison between our measurements and 



the results of quantum chemical calculations allowed us 
to characterize the lowest-energy transitions. Based on our 
assignments, we constructed an exciton coupling model to 
describe the absorption spectrum of 6-MI dinucleotide in 
B-form double-stranded (ds) DNA. The results are in 
good agreement with experimental data. 

MATERIALS AND METHODS 

Experimental procedures 

Chemicals 

6-methylisoxanthopterin riboside (r(6-MI)) was obtained 
from Fidelity Systems, Inc. and used without further puri- 
fication. PVA (M w = 85 000-124000) was purchased from 
Sigma Aldrich and used as received. All organic solvents 
were spectrophotometric grade. Solution phase absorp- 
tion spectra were measured using 13.2 uM r(6-MI) 
dissolved in 5mM aqueous sodium phosphate buffer, 
pH 7. 

Film preparation 

A 9.1% (by weight) aqueous PVA solution was prepared 
by dissolving the polymer in nanopure water while heating 
to 95°C, and then allowing the solution to cool to room 
temperature. Sample solutions were prepared by 
combining 6 ml of the aqueous PVA solution with 
1.33 ml of 0.45 mM r(6-MI) in methanol, and stirring the 
mixture for several minutes. A reference solution was 
prepared by following an identical procedure but with- 
holding the r(6-MI). The sample and reference solutions 
were then placed under vacuum for 1 h in a desiccator that 
contained a small amount of water. The solutions were 
transferred to Eppendorf tubes and spun down in a cen- 
trifuge at 10 000 r. p.m. for 15 min to remove bubbles and 
non-dissolved PVA. Films were cast by pouring 3 ml of 
sample (or reference) solution into 2.5 x 2.5 x 0.5 cm 3 
wells, which were cut into paraffin wax molds supported 
on a glass sheet. The films were then allowed to dry for 2 
days in an area free from dust and air currents. Films were 
stretched to 1.5x their original length while they were 
heated using a hot-air gun. 

Spectroscopy 

Solution phase absorption spectra were measured using 
a Cary 3E UV-visible spectrophotometer. LD and 
absorption spectra of PVA films were simultaneously 
measured using a Jasco CD spectrophotometer running 
in LD mode (Horiba Jobin-Yvon, J-720). Sample LD 
spectra were corrected for background from the PVA 
host material by subtracting the LD of a stretched refer- 
ence film from that of the stretched sample films. 
Significant host absorption was observed in sample films 
at wavelengths below 250 nm. The reduced linear dichro- 
ism (LD r ) spectra were calculated using the absorption 
spectrum of r(6-MI) in aqueous buffer, which was scaled 
so that the intensity of its lowest-energy absorption peak 
(where PVA absorption is negligible) matched that of 
r(6-MI) in the film. 
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Quantum chemical calculations for electronic excitation 
energies and electric dipole transition moments 

Low-lying electronic excitation energies and EDTMs of 
small organic molecules can be predicted using modern 
quantum computational methods with reasonable accur- 
acy (25-29). The computed excitation energies are typic- 
ally accurate to within 0.2-0.4 eV (ca. 1600-3200 cm" 1 ) for 
low-lying valence electronic excitations (25,26,29). The 
accuracy of computed EDTMs and oscillator strengths 
is less studied (28). However, the shapes of electronic ab- 
sorption spectra predicted by quantum chemical methods 
tend to match experimental absorption spectra quite well, 
especially in their low-energy parts. 

In this work, we use three well-established quantum 
chemical methods to perform ground-state geometry op- 
timizations and to compute excitation energies and tran- 
sition moments for the 10 lowest electronic excitations of 
N8-methyl-6-MI, namely approximate second-order 
coupled-cluster with resolution of the identity approxima- 
tion (RICC2) (30,31) and density functional theory using 
the Becke 3-parameter hybrid (B3LYP) (32) and the <bB97 
long-range corrected functional (33). Triple-zeta valence 
basis sets with two sets of polarization functions 
(TZVPP) (34) and corresponding auxiliary basis sets (35) 
were used throughout. The stabilities of the optimized 
ground-state structures were established by force 
constant calculations. Owing to pyramidalization of the 
NH 2 group, the ground-state structure of N8-methyl-6- 
MI belongs to the CI point group. RICC2 and B3LYP 
calculations were performed using the TURBOMOLE. 
program package (36). coB97 calculations used the 
Q-Chem program (37). 

Moments of inertia of r(6-MI) were computed using 
B3LYP-optimized ground-state structures (32). The 
resolution-of-the-identity approximation was used for 
the Coulomb energy (RI-J) (38). TZVPP basis sets were 
used along with corresponding auxiliary basis sets. To de- 
termine the lowest-energy conformation of r(6-MI), differ- 
ent conformations of the ribose ring and the glycosidic 
substituent were sampled by the Dreiding force field (39) 
and subsequently optimized using the B3LYP functional. 

Theoretical considerations for linear dichroism 
measurements 

LD experiments were performed on r(6-MI), in which the 
ribose sugar is covalently bound at the N8 position of the 
6-MI free base (see Figure 1). Samples were prepared by 
suspending the nucleoside in PVA films. By applying a 
uniaxial strain to these films, a net alignment of the exci- 
tation transition moments was established (19,23,40,41). 
We performed LD measurements on both stretched and 
un-stretched PVA films. 

LD is the difference between absorption of light 
polarized parallel and perpendicular to an orientation 
axis, which is defined as the direction of the uniaxial 
strain applied to the PVA films. 

LD(y) = ^||(v) - A ± (v) (1) 



The absorption of each polarized component is propor- 
tional to the square of the projection of the transition 
moment vector onto the electric field vector. The extents 
to which the parallel and perpendicular transition moment 
projections differ — and hence, the magnitude of the LD 
signal — depends on both the orientation of the r(6-MI) 
transition moments and on how narrowly distributed the 
r(6-MI) orientations are in the film. The reduced linear 
dichroism (LD r ) is defined as 

LD r (y) = LD(v)/A iso (v) (2) 

where A iso (v) = l\A\i(v)+2Ax(v)\ is the corresponding ab- 
sorption of light in an isotropic sample. The quantity 
LD r (y) depends only on the orientations of the transition 
moments that absorb at frequency v, and not on their 
strengths. For a purely polarized transition, the value of 
LD r is expected to be constant across the absorption 
band. 

The distribution of r(6-MI) orientations in the film is 
specified by the Saupe orientation parameters S xx , S yy and 
S zz , where the z-axis is taken to be the uniaxial direction of 
the film (19). For a planar molecule with in-plane and 
out-of-plane transition polarizations, the x-axis is 
assumed to be perpendicular and the y- and z-axes are 
parallel to the molecular plane. The relationship between 
the value of LD] and the angle of the ith transition 
moment vector relative to the film orientation axis is 
given by 

LD- = 3(S vv sin 2 #,+S z _cos 2 0,) (in-plane transitions) (3) 

Ideally, the Saupe parameters would be determined in 
an independent experiment. For example, the infrared LD 
spectrum could be used, in principle, to solve Equation (3) 
for the Saupe parameters under the assumption that the 
molecule aligns with its longest axis parallel to the film 
stretch direction, and that the vibrational transition 
moments of the carbonyl groups can be clearly identified 
(22,24). Unfortunately, our attempts to measure the infra- 
red LD spectrum of our samples were unsuccessful 
because the limited solubility of r(6-MI) resulted in insuf- 
ficient signal intensity compared with background absorp- 
tion by the PVA host material. We therefore followed an 
alternative approach to estimate the Saupe parameters by 
using the following relations, where S xx is the orientation 
parameter for the out-of-plane polarized x-axis (19,24). 

LD- = 3S XX (out-of-plane transitions) (4) 

3S yy < LD" < 3S :: (in-plane transitions) (5) 

and 

S XX + Syy + S ZZ = 0 (6) 

In the approach we followed, the value of S xx was 
estimated from the average of the LD r values correspond- 
ing to the out-of-plane polarized transitions. The values of 
Syy and S zz were estimated from the minimum and 
maximum LD r values, respectively, of in-plane polarized 
transitions. Equation (6) was used to calculate S :: because 
a transition with positive LD r was not available from 
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Calculations: R : 
Experiments: R : 

Figure 1. The structure of the fluorescent guanine analogue 6-MI is 
shown with the atomic numbering system used in this work. The 
in-plane angle 5 is defined relative to the C9-C10 axis, as shown. For 
our LD measurements on the 6-MI nucleoside (r(6-MI)), the R-group 
was ribose. Quantum chemical calculations were carried out on the 
8-methyl-substituted 6-MI. 



which to estimate it. The Saupe parameters were then used 
in Equation (3) to calculate the value of for the /th 
in-plane polarized electronic transition, which yielded 
two possible values of 8 h one positive and one negative. 
We selected between these possibilities by comparison 
with the computational results. Finally, we used the 
computed ground-state geometry of r(6-MI) to calculate 
its minimum principal moment of inertia axis, which lies 
approximately along the longest axis of the molecule. This 
was assumed to lie parallel to the orientation axis of the 
stretched films. The angle 0 t was then related to the angle 
Si, which describes the orientation of the EDTM of the /th 
transition polarized in the plane of the molecule. 



beneath the shoulder of the feature centred at 
47 547 cm -1 (peak 6). We determined the values of LD r 
for each transition using the method of Thulstrup, Eggers, 
and Michl (TEM method) (42), which involves making 
linear combinations of the polarized spectra, 
A\\(y) — dA±(v) for different values of the reduction coef- 
ficient d. The value dj for which the /th spectral feature 
vanishes yields the corresponding value LD. of the tran- 
sition, according to 



, dj-\ 
' d,+2 



(7) 



We performed a global numerical fitting analysis of the 
experimental absorption spectrum and the LD r spectrum, 
in which each of the six absorptive features was modelled 
as a Gaussian line shape with variable frequency, line 
width and intensity. This was accomplished by minimizing 
the target function described by Equation (8) with the aid 
of the nonlinear global optimization package KNITRO 
(43), where A ca ] c (v) is given by Equation (9) and 
LD r calc (v) is given by Equation (10). 

50 000 cm- 1 

x 1 = ^ex P (?) - A cak {v)] 2 +[LD^ p {v) - LDr calc {v)f 



i=25 000 cm- 



6 (v-5,-) 2 



A C alc(v) = ^2 he 2 °: 



(8) 
(9) 



RESULTS 

The absorption spectrum of r(6-MI) in aqueous buffer 
exhibited three major peaks centred at 29 687, 34 596 
and 47 547 cm" 1 (337, 289 and 210 nm), which we 
identify as peaks 1, 2 and 6, respectively (see Figure 2A). 
The absorption spectrum of r(6-MI) in PVA films ex- 
hibited the same features as in isotropic solution, 
although the highest-energy transition was obscured by 
background absorption of the PVA host material. The 
LD spectrum of r(6-MI) in un-stretched PVA films was 
weak. In stretched PVA films, the LD spectrum exhibited 
features at the same energies as in isotropic solution. 
However, as expected, the relative intensities of each ab- 
sorptive feature in the stretched films deviated from those 
measured in isotropic solution. Four films were examined, 
and the analysis below was carried out using the average 
of two LD spectra of one representative film. 

We constructed the experimental reduced linear dichro- 
ism (LD 1 ) by dividing the LD obtained from the stretched 
PVA films by the absorption measured in aqueous buffer, 
according to Equation (2) (see Figure 2B). We chose to use 
the buffer solution spectrum rather than that of the 
un-stretched PVA films in Equation (2) because of the 
background absorption of the PVA at wavelengths 
shorter than 250 nm. On examination of the LD r , it was 
revealed that three additional transitions occur over the 
range 38 000^44 000 cm" 1 (labeled as peaks 3, 4 and 5). 
These features were initially not evident from the experi- 
mental absorption spectrum because they were buried 



(10) 



1=1 



This procedure yielded optimized parameters for each 
of the six transitions. We note that our results were inde- 
pendent of film preparation conditions (using methanol or 
dimethyl sulfoxide (DMSO) as the solvent for r(6-MI)) 
and independent of whether the film or solution absorp- 
tion spectrum was used to calculate the LD r spectrum. 
The results of the fitting analysis are shown in Figure 
2 A and are listed in Table 1. 

We next compared the results of the quantum chemical 
calculations with those of the experimental absorption and 
LD r spectra (summarized in Table 2). We noted that there 
are significant differences between the LD r values of peaks 
1, 2 and 6 and those of peaks 3, 4 and 5, despite the ab- 
sorbance being weak in the region of the latter group, 
suggesting the presence of hidden weak transitions with 
significantly different orientations from the in-plane 
polarized transitions. The RICC2 results indicated the fol- 
lowing progression in order of increasing energy; two 
strong in-plane polarized transitions, which we assigned 
to peaks 1 and 2 of the experimental spectra; four weak 
out-of-plane polarized transitions, of which only peaks 3, 
4 and 5 could be identified in the data and a series of 
strong in-plane polarized transitions, of which experimen- 
tal peak 6 was assigned to the lowest in energy. The data 
therefore matched the ordering of transitions found with 
the RICC2 method, provided that peaks 3, 4 and 5 were 
assigned to out-of-plane polarized transitions. 
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Figure 2. (A) The experimental absorption spectrum of r(6-MI) in aqueous buffer (black) is shown in comparison with a fit (grey) of the data to the 
sum of six Gaussian functions (dashed grey). Transitions are numbered according to the same scheme as in Table 1. (B) The experimental LD r 
spectrum of r(6-MI) in the stretched PVA film (black) is shown in comparison with the calculated LD r spectrum (grey), which is based on the 
experimental transition moment orientations and intensities, energies and peak widths obtained from a simultaneous fit to the data shown in (A) and 
(B). (C) Plot (A), shown as a function of wavelength in nanometers. (D) Plot (B), shown as a function of wavelength in nanometers. 



Table 1. Results of global fitting of r(6-MI) absorption and LD r 
spectra 



Transition 


s (NT 1 cnr 1 )' 1 


v (cirr') b 


X (nm) c 


a (cm -1 ) d 


1 


14400 


29 687 


336.8 


1815 


2 


11200 


34 596 


289.1 


1452 


3 


2900 


38 327 


260.9 


1344 


4 


3300 


41282 


242.2 


1371 


5 


8300 


43 923 


227.7 


2107 


6 


32 000 


47 547 


210.3 


1952 



''Molar extinction coefficient. 
b Centre transition energy in wave numbers. 
c Centre transition wavelength in nanometers. 
d Absorption line width. 



RICC2 predicts electronic excitation energies that are 
ca. 0.35 eV (ca. 2800cm -1 ) too high for the strong peaks 
1 and 2, consistent with previous benchmark studies (27). 
The accuracy of theoretical predictions for higher-lying 
electronic excited states is less certain. The excitation 
energy for peak 6 is underestimated by ca. 0.3 eV (ca. 
2300 cm -1 ) by RICC2 compared with experiment. The 
coB97 functional systematically overestimates the experi- 
mental excitation energies for the peaks 1-5 by ca. 0.81 eV 
(ca. 6800 cm -1 ), which is similar to the 0.5^3.6 eV error 
found in a recent benchmark study (28). The computed 



oscillator strengths of peaks 1 and 2 are in good agreement 
with experiment for both RICC2 and coB97 methods. In 
comparison, the oscillator strength of peak 6 is consider- 
ably underestimated by RICC2 and is slightly 
underestimated by coB97. These results suggest that peak 
6 is likely to be prone to mixing with other high-lying 
excited states. As a result, the predicted excitation 
energies and transition moments are expected to be less 
accurate. 

We estimated the value of S xx to be —0.157 according to 
Equation (4) and by using the average of the LD r values 
corresponding to peaks 3, 4 and 5. Of the remaining 
in-plane polarized transitions, peak 1 had the lowest 
value of LD r (—0.268), which was used in Equation (5) 
to place an upper limit on the value of S rr (<— 0.089). 
Because no transition was observed with a positive value 
of LD 1 , we estimated a lower limit on the value of S zz 
(>0.246) using Equation (6). A lower limit on the value 
of Syy was estimated using the properties of an orientation 
triangle (see Figure 3) (19). The orientation triangle de- 
scribes the relative magnitudes of the Saupe orientation 
parameters for molecules of different symmetries. 
Molecules with ideal rod-like symmetry are expected to 
fall along the lower side of the triangle, whereas molecules 
with ideal disk-like symmetry fall along the upper side of 
the triangle. Molecules of intermediate symmetries must 
fall within the boundaries of the triangle. Using the value 
we previously determined for Sxx = —0.157, we note that 
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Table 2. Computed and experimental electronic excitation parameters for 6-MI 



Tran . 




RICC2 






coB97 






Experimental 




v (cm ') a 


/ 


5(°) c 


v (cm ') 


/ 


5 (°) 


v (cm ) 


f d 


5, (°) e 


5 2 ("/ 


r 


32 342 


0.031 


76 


35 659 


0.001 


oop 










i 


32410 


0.28 


76 


36372 


0.35 


73 


29 687 


0.28 


66 


8S 


2 


37 756 


0.12 


-22 


41287 


0.15 


-28 


34 596 


0.18 


-61 


-69 


3 


40 115 


0.0001 


oop 


45 028 


0.0005 


oop 


38 327 


0.042 


oop 


oop 


4 


42 372 


0.0002 


oop 


45 981 


0.0001 


oop 


41282 


0.049 


oop 


oop 


5 


42 592 


0.0003 


oop 


50 796 


0.0019 


oop 


43 923 


0.19 


oop 


oop 


5' 


42 994 


0.001 


oop 
















6 


45 158 


0.16 


-49 


52816 


0.51 


-61 


47 547 


0.68 


-66 


-78 



"Centre transition energy in wave numbers. 
b Oscillator strength. 

"Angle of in-plane polarized EDTM relative to C9-C10 axis in degrees (oop = out-of-plane). 
d Oscillator strength obtained by integrating the Gaussian functions given in Table 1. 
e EDTM angle corresponding to orientation parameters obtained from orientation triangle. 
f EDTM angle corresponding to orientation parameters obtained from Equation (5). 




Figure 3. An orientation triangle for anisotropic molecules in stretched 
polymer films. The large dot specifies the in-plane orientation param- 
eters obtained for 2-aminopurine nucleoside [24]. The upper small dot 
indicates the orientation parameters corresponding to the largest 
possible value of S yy that is consistent with the experimental LD r 
data for r(6-MI). Because S xx ( = —0.157) is defined by the experimental 
LD r s of the out-of-plane polarized transitions, the values of S yy and S- z 
must fall on the line: S y) ,+S z , = —S xx [see Equation (6)]. The intersec- 
tion of this line (dashed grey) with the boundary of the orientation 
triangle (lower small dot) yields the smallest value of Syy, and the 
largest value of SU that are consistent with the experimental data. 

the values of S rr and S zz must fall on the line 
Syy+Szz = —S xx , as required to satisfy Equation (6). The 
intersection between this line and the boundary of 
the orientation triangle provides the lowest possible 
value of Syy ( = —0.157) and the highest possible value 
of S zz ( = 0.314) that are consistent with the experimental 
data. The absolute value of 6, was then calculated for 
each transition using Equation (3). The order param- 
eters and values of 9j for in-plane polarized transitions 1, 
2 and 6 obtained by this analysis are summarized in 
Table 3. 



Table 3. Experimental limits on orientation parameters and in-plane 



polarized EDTM an 


lies 












s xx 


Syy S 


\0i\ (°) 


m (°) 


1% 1 (°) 


Orientation triangle 


-0.157 


-0.157 0.314 


67.7 


59.8 


64.5 


Equation (5) 


-0.157 


-0.089 0.246 


90.0 


67.0 


75.9 



We determined the minimum principal moment of 
inertia axis from B3LYP calculations of the optimized 
ground-state geometry of r(6-MI), finding that it lies 
roughly parallel to the glycosidic bond. By assuming 
that r(6-MI) orients with this axis parallel to the PVA 
film stretching direction, we converted the lab-frame 
angles to the angles S h which describe the orientations 
of in-plane polarized EDTMs relative to the C9-C10 axis 
(the two central carbons of the 6-MI base, see Figure 1). 
For each electronic transition, we obtained upper and 
lower bounds to the angles Sj by estimating the range of 
orientation parameters, as described above (see Figure 4). 
The RICC2 prediction for the lowest-energy transition 
(peak 1) falls within the experimental limits, but this is 
not the case for the next two higher-energy in-plane 
polarized transitions (peaks 2 and 6). As discussed 
above, state mixing affects the computed orientations of 
higher-energy transition moments more than low-energy 
ones. However, we note that the prediction to peak 6 by 
the coB97 method falls only 5° outside of the experimen- 
tally determined range, and the ordering of the transitions 
from 'most j-polarized' to 'most z-polarized' are in agree- 
ment between experiment and both RICC2 and coB97 
methods. Our experimental and computed results for the 
orientation of the EDTM of peak 1 are consistent with the 
orientation computed using the B3LYP functional and 
6-311 +G(d,p) basis set in reference 8 (8 = 79°) (8). 

DISCUSSION 

As mentioned in the Introduction to this article, 6-MI can be 
used as a fluorescent guanine analogue probe to study the 
local conformations of nucleic acids. A number of 
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Transition 1 




Experimental range 

RICC2 

coB97 

B3LYP 

Figure 4. The results of experimental and calculated in-plane polarized 
EDTMs of 6-MI are shown. Grey hourglass regions indicate the range 
of EDTM orientations corresponding to the experimental limits for S yy 
(indicated by upper and lower small dots in Figure 3). The green 
(RICC2), dark blue (coB97) and light blue (B3LYP) arrows indicate 
the EDTMs obtained from calculations. (A) Transition 1; 
(B) Transition 2; (C) Transition 6. 

spectroscopic observables can yield qualitative information 
about local base environment, without detailed knowledge 
of the EDTMs in 6-MI. For example, 6-MI fluorescence 
intensity can indicate the extent of local base stacking in 
the nucleic acid or protein-nucleic acid complex being 
studied. Spectroscopic measurements that are sensitive to 
the exciton coupling between the adjacent 6-MI residues, 
such as CD and 2D FS, have the potential to provide 
more quantitative information about local base conform- 
ation, if we have accurate knowledge of the EDTMs. 

In 6-MI dinucleotide-substituted DNA constructs, 
exciton coupling is sensitive to the magnitude, separation 
and relative orientation of the 6-MI EDTMs. In the 
point-dipole approximation (16), the exciton coupling 
between two 6-MI chromophores can be described using 
Equation (1 1), where fi„ is the transition dipole of the «th 
[n e {1,2}] chromophore, R is the inter-chromophore 
centre-to-centre vector and So is the vacuum permittivity. 

F|2 = 4^ /i '- (1 - 3 ^ ) - M2 (U) 



Thus, a determination of the coupling for specific con- 
formations requires knowledge of the EDTM orientation 
within the 6-MI chromophore. 

To test our determination of the lowest-energy EDTM 
of 6-MI, we modelled the exciton coupling between these 
transitions in a 6-MI di-probe-substituted dsDNA con- 
struct (3) (see Figure 5). The dipole strength |/x| 2 of the 
transition can be determined from the absorption 
spectrum of a mono-probe-substituted dsDNA construct, 
according to Equation (12). 

oc 

\fi 2 =— — - / dv— (12) 

— 00 

Here, a is the absorption coefficient, ti is Planck's 
constant divided by 2ji, c is the speed of light and N A is 
Avogadro's number. The factor f^Ldv^ is the 
integrated optical line shape of the lowest-energy transi- 
tion, measured in wave numbers and divided by its fre- 
quency. We fit the lowest-energy absorption peak of a 
6-MI mono-probe-substituted dsDNA construct to a 
Gaussian line shape, as shown in Figure 5A. We thus 
determined the line width, intensity and transition 
energy of the 6-MI monomer probe within the dsDNA 
construct to be ofwhm = 3480cm -1 , 7=0.011 absorb- 
ance units/uM 6-MI and Ei = 28 183 cm -1 , respectively. 
Using Equation (12), we estimated the norm of the tran- 
sition moment of the 6-MI monomer to be |/x.| = 3.65 D. 
This is to be compared with a computed value of 4.27 D, 
which was obtained using the RICC2 method. 

Using the values determined experimentally for the 
monomer transition energy and dipole strength, we 
calculated the absorption spectrum expected for the 
6-MI dinucleotide-substituted B-form dsDNA construct 
(3) (see Figure 5B). These calculations assumed the 
in-plane orientation of the EDTM that we calculated 
using the RICC2 method (8 = 76°), which falls within 
the boundaries established by our experiments. We 
further assumed that the coordinates of the two adjacent 
6-MIs were the same as those of two guanines in B-form 
dsDNA, which we obtained from crystallographic data 
(protein data bank reference code 1BD1). 

Based on the exciton coupling model described above 
for the 6-MI dinucleotide dsDNA construct, we obtained 
values for the coupling V i2 = 814 cm -1 and relative dipole 
angle a = 36.9°. The absorption spectrum of the di-probe 
DNA construct was calculated assuming that the exciton 
coupling between 6-MI residues leads to a splitting of the 
transition energies with e± = Si ± Vn and intensities 
7± = (1 ± cosa). Because the adjacent bases in 
duplex DNA adopt a face-to-face geometry, the sign of 
V\2 is positive, such that the resulting spectrum exhibits a 
strong, blue-shifted transition and a weaker red-shifted 
transition, as shown in Figure 5B. We note that the agree- 
ment between the calculated and experimental spectra is 
excellent. For different possible values of the 
lowest-energy EDTM angle S, the magnitude of the 
exciton coupling varies from 136 to 856 cm -1 , and only 
values of 8 near our experimentally determined and 
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6-MI mononucleotide substituted dsDNA 
A 5' CTAATCATTGTgGTTCGGTCCTTGC C 
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B 6-MI dinucleotide substituted dsDNA 

5' CTAATCATTGTggTTCGGTCCTTGC 
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Figure 5. (A) A fit is shown of the low-energy region of the absorption spectrum of dsDNA containing one 6-MI base substitution (black) to a 
Gaussian (purple). The fit was optimized using only the 23 310-29 326 cm -1 region of the absorption spectrum to model only the lowest-energy 6-MI 
transition (transition 1 in our assignments), and to avoid effects from overlapping transitions in the native DNA bases. (B) The determined strength 
and orientation of the lowest-energy EDTM in 6-MI was used to calculate the coupling between adjacent 6-MIs expected for a B-form dsDNA 
structure. The spectrum of dsDNA containing two adjacent 6-MI residues (black) is shown along with the spectrum modelled using the calculated 
coupling and crystallographic structural parameters (purple, underlying exciton transitions in dashed blue and dashed red). The sequence of each 
6-MI containing DNA strand is shown above its spectrum. The complimentary strand is not shown. (C) Plot (A), shown as a function of wavelength 
in nanometers. (D) Plot (B), shown as a function of wavelength in nanometers. Experimental data from reference 3. 



computed values reproduce the observed shift in the ab- 
sorption spectrum on addition of a second 6-MI residue. 

While the calculated spectrum shown in Figure 5B 
captures the observed blue shift of the absorption 
maximum, there is a significant low-energy shoulder in 
the experimental spectrum that is not reproduced by the 
exciton coupling model. The calculated spectrum has only 
a slight shoulder in this region owing to the presence of a 
weak red-shifted exciton transition. We note the presence 
of a similar shoulder in the absorption spectrum of the 
6-MI mono-probe-substituted DNA construct, which rep- 
resents a deviation of the monomer line shape from 
Gaussian. A similar shoulder is visible in the 
low-temperature absorption spectrum of r(6-MI) and is 
likely due to an underlying vibronic progression (8). 



CONCLUSIONS 

Using a combination of absorption and LD measurements 
of 6-MI in stretched PVA films, we identified six electronic 
transitions that occur within the spectral region between 
25 000 and 50 000 cm -1 . The range of possible orientations 
for each in-plane transition moment was determined ex- 
perimentally. We performed quantum chemical calcula- 
tions (RICC2 and coB97) to predict (in order of 



increasing energy) two in-plane transitions, three weak 
out-of-plane polarized transitions and a strong in-plane 
polarized transition. We found that the orientation of 
the lowest-energy transition moment, which was predicted 
by both RICC2 and coB97, agreed well with experiment. 
We compared the absorption spectra of two 6-MI- 
substituted dsDNA constructs, each labelled with just 
one or two adjacent 6-MI residues. We found that the 
strength and orientation of the lowest-energy EDTM 
determined in this work accurately predicts the experimen- 
tally observed blue shift in the absorption due to exciton 
coupling on addition of a second 6-MI residue. In future 
work, these results will be used to interpret in structural 
terms 2D FS spectra reflecting the conformations of 
site-specifically placed 6-MI probes in biologically active 
protein-nucleic acid complexes. For example, such studies 
have the potential to reveal valuable insights into the role 
of local base pair conformation and breathing fluctuations 
on the biochemical pathways that govern DNA replica- 
tion and transcription. 
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